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Heavy quark potential in the instanton liquid model 
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We study the heavy quark potential in the instanton liquid model by carefully measuring Wilson loops out 
to a distance of order 3/m. A random instanton ensemble with a fixed radius p — 1/3/m generates a potential 
V(R) growing very slowly at large R. In contrast, a more realistic size distribution growing as p 6 at small p and 
decaying as p~ 5 at large p, leads to a potential which grows linearly with R. The string tension, however, is only 
about 1/10 of the phenomenological value. 
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1. INTRODUCTION 

There is growing evidence that instantons play 
a very important role in QCD. The development 
of the Interacting Instanton Liquid Model Q al- 
lows one to calculate nonperturbatively to all or- 
ders in the 't Hooft interaction, and results show 
that it correctly generates the quark condensate 
and most salient properties of light hadrons. As 
reviewed at this conference 0, the essential fea- 
tures and parameters of this model have now been 
confirmed on the lattice, and there is strong lat- 
tice evidence that instantons and their associated 
zero modes play a significant role in hadron struc- 
ture. 

Motivated by a provocative result by Fuku- 
shima et al. J3| reporting a string tension in the 
instanton liquid model close to the phenomeno- 
logical value, we have undertaken a careful calcu- 
lation of the instanton induced heavy quark po- 
tential. 



2. SIMULATION DETAILS 

The gauge field of a single instanton in SU(2) 
in the singular gauge (centered at the origin) is 
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given by 



Afi{x) — 2T a r) a ^ v — 



x 1 x 2 + p z 



(1) 



where f) apiV is the 't Hooft symbol. Replacement 
of fj af j,u by rj aflL/ yields the solution for an anti- 
instanton. The multi-instanton gauge configura- 
tion in the dilute instanton liquid model is gener- 
ated by a linear superposition of the fields from 
a distribution of instantons and anti-instantons. 
Each instanton and anti-instanton has a random 
center location and is embedded in SU(3) with a 
random color orientation. 

We consider two distributions of the instanton 
size, p, for 517(3): (1) a fixed size distribution, 
with p = 1/3 /to; and (2) a phenomenological 
distribution, with D(p) ~ p 6 at small p as ex- 
pected from the dilute instanton gas approxima- 
tion and D(p) ~ p~° at large p as would occur if 
the perturbative running of the coupling constant 
g were frozen at sufficiently large p. We use the 
parameterization by ShuryakQ] for this variable 
size distribution. 



D(p) = C- 



0o - 5 + p 3 - 5 ) — 



(2) 



We choose po so that the average instanton size 
p = 1/3/m and C is a normalization constant. 

We generate ensembles of multi-instanton con- 
figurations in an open box of size (L x , L y , L z , L t ) 
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Figure 1. V(T,R) as a function of T for the 
fixed instanton size distribution p = 1/3/m at 
an instanton number density n — 1.0/m~ 4 . Each 
data set from bottom to top corresponds to R = 
0.4/m, 0.8/to, 1.2/to, 1.6/m and 2.0/rn, respec- 
tively. 



Figure 2. V(T, R) as a function of T for a variable 
instanton size distribution described in Eq. (0) 
at an instanton number density n — l.Ofm . 
Each data set from bottom to top corresponds to 
R = 0.4/m, 0.8/m, 1.2/m, 1.6/m and 2.0/m, 
respectively. 



(13.6/to, 7.2/m, 7.2/m, 20/m). Wilson loops 



W(-R, T) = trP exp[i i A^dxJ 

JRxT 



(3) 



are measured in a 6.4/to x 12.8/m rectangle in 
the central plane. This leaves a distance of at 
least 3.6/m or more than 10 times the average 
instanton size p from the measured Wilson loop 
to the 4-d boundary so that the edge effects from 
the open boundary are small. The largest Wilson 
loop is 3.2/m x 6.4/m in size. The discretization 
Aa: in the path ordered integral is 0.05/m, which 
is about 1/6 of the average instanton size. 

We study three different instanton number den- 
sities, n = N/V = 0.5/m~ 4 , 1.0/m~ 4 and 
1.5/m~ 4 for both fixed and variable instanton 
size distributions. The number of gauge configu- 
rations studied for each different n is 1600, 8000 
and 1600, respectively. 



3. RESULTS 

The fundamental difficulty in this calculation 
is measuring with adequate statistical accuracy 
the large Wilson loops required to determine the 
potential at large distances. We determine the 
heavy quark potential V(R) from rectangular 
loops of different time extent as follows: 



V(R) = lim V(T,R). 



(4) 
(5) 



To obtain the most accurate measurement of the 
potential, we plot V(T, R) vs. T and find the 
plateau corresponding to V(R). 

Figures |l| and || show such plots at different 
R for both fixed and variable instanton size dis- 
tributions at an instanton number density n = 
1.0/m -4 . For the fixed instanton size distribu- 
tion, p = 1/3 fm, plateaus for V(T,R) set in at 
about T > (1 ~ 2)R. We have good plateaus for 
all R < 3.2/to at n = 1.0/to~ 4 . In contrast, for 
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Figure 3. Heavy quark potential in the instanton 
liquid model for a fixed instanton size p = 1/3/m. 
Three different instanton number densities are 
plotted. 



Figure 4. Heavy quark potential for a variable in- 
stanton size distribution with p — 1/3 fm. Three 
different instanton number densities are plotted. 



the variable instanton size distribution, not only 
do the plateaus set in at a larger T > (2— 3)R, but 
also the statistical fluctuations are much larger 
at large R, especially for R > 2.0/m. ft is clear 
that this behavior of the variable size distribu- 
tion arises from large instantons in the tail of the 
distribution. Even with an ensemble as large as 
8000 configurations, we are only able to obtain a 
potential V(R) for R < 2.0/m at n = 1 .0/m~ 4 . 

Figures [| and |I] show the heavy quark poten- 
tial out to the maximum distance we can reliably 
measure. The salient results are the following. 
At short distance, the slope is identical for both 
the fixed and the variable size distributions at a 
given n and is proportional to n. At the value 
n = 1.0 fm of the instanton liquid model, this 
slope is ~ O.f GeV/ fm, corresponding to 1/10 of 
the physical string tension. We note that by scal- 
ing, the physical string tension would be obtained 
at n = 1.0/m -4 by increasing the mean value of p 
by f 1 / 4 to 0.59 fm. When R is much greater than 
the largest instantons in the medium, we expect 
V(R) to approach a constant [Q . For the distribu- 
tion in Eq. (S), however, it is interesting that the 



potential is essentially linear in the region of 0/m 
- 2/m. Hence, even though the random instan- 
ton liquid is not strictly confining, in the region 
of physical interest for hadron structure, it has a 
significant linear component. 
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